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We show that a flow (timelike congruence) in any type B\ warped product spacetime is uniquely 
and algorithmically determined by the condition of zero flux. (Though restricted, these spaces 
include many cases of interest.) The flow is written out explicitly for canonical representations of the 
spacetimes. With the flow determined, we explore an inverse approach to Einstein's equations where 
a phenomenological fluid interpretation of a spacetime follows directly from the metric irrespective 
of the choice of coordinates. This approach is pursued for fluids with anisotropic pressure and shear 
viscosity. In certain degenerate cases this interpretation is shown to be generically not unique. The 
framework developed allows the study of exact solutions in any frame without transformations. We 
provide a number of examples, in various coordinates, including spacetimes with and without unique 
interpretations. The results and algorithmic procedure developed are implemented as a computer 
algebra program called GRSource. 
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I. INTRODUCTION 



The usual procedure for finding exact solutions of Ein- 
stein's equations involves writing down a phenomenolog- 
ical energy-momentum tensor, often a perfect fluid, in a 
set of coordinates, frequently comoving, so that the field 
equations can be integrated, often with the aide of simpli- 
fying assumptions Q ■ In view of the difficulty of solv- 
ing Einstein's equations, inverse problems are important 
[Sj . An inverse problem of interest (not involving the clas- 
sification problem for the Ricci tensor explicitly) may 
be stated thus: Given a spacetime {M., g) with manifold 
M. and Lorcntzian metric g, what, if any, fluid flow could 
generate {M. , g) via Einstein's equations? This question 
is explored here for the case of type Bi warped product 
spacetimes (see below). The problem is of interest since, 
as is shown, the flow (velocity field) is uniquely deter- 
mined subject to a zero flux condition. With the flow 
determined, physical parameters, subject to a fluid de- 
composition that includes anisotropic pressure and shear 
viscosity, can be extracted directly from the metric ir- 
respective of coordinates. The framework developed in 
this paper allows the study of exact solutions in any 
frame without transformations. Surprisingly little work 
has been done in non-comoving frames |^ . The invariant 
procedure developed is algorithmic and suited to com- 
puter algebra projects such as GRDB ^ where, given a 
metric, one might like to know if it is necessarily, for ex- 
ample, a perfect fluid, or on the contrary, incompatible 
with a perfect fluid. The paper is organized as follows: 
In section II, we explore the zero flux condition and in 
section III we derive explicit forms of the velocity field for 
canonical coordinate types. In section IV we explore the 



phenomenology of the non-conducting fluid source show- 
ing that degenerate cases exist. Examples that illustrate 
the power of the results obtained are given in section V. 
Section VI is a summary. 



II. ZERO FLUX 

We consider warped product spacetimes of class Bi 
0,0. These can be written in the form 



ds 



M 



ds'^^ {x' , x') + C{x°')dsi^^ (a;^ x*) 



(1) 



where C{x°') = r{x^ , x^yw{x'^ , x'^y , sigCSi) — and 
sig{'S2) = 2e (e = ±1)- Although very special, these 
spaces include many of interest, for example, all spheri- 
cal, plane, and hyperbolic spacetimes. We write 



ds|^ = a{dx'^f + 2bdx^dx^ + c{dx^f 



(2) 



with a,b and c functions of (a;^,a:^) only. A congruence 
of unit timelike vectors (a "flow" in what follows) = 
(u^,u'^,0,0) have an associated unit normal field n" (in 
the tangent space of Ei) satisfying Uau" = O,^^^" — 1 
[3. It follows that Ua = i){x^,x'^){u'^,-u^,0,0) {ip a 
normalization factor). The timelike condition on is 

-l^a{u^f + 2bu^u^ + c{u^f (3) 

and {Ai , g) is time orientated by the restriction 

> 0. (4) 

The condition 

G^^M^n^ = 0, (5) 
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where is the Einstein tensor, can be written as 

Au^u^ - B{u^)^ + C{uY ^0, (6) 
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where 

A = G\-Gl B = Gl C = G\. (7) 

Condition ((SJ is the defining criterion for the flows con- 
sidered here. In general, A, B and C are functions of 
through a;'*, assumed not to vanish simultaneously. 
Equations Q and © determine the flow uniquely. 
The flow need not exist. For example, \i A = G = and 
B 7^ then there exists no such flow. The assumption of 
"comoving coordinates" {v? = in the present notation) 
imposes restrictions on the coordinatization of Si, but it 
is already clear from ^ that the existence of such a flow 
requires B = Q. 



III. u° FOR CANONICAL REPRESENTATIONS 

OF El 

There are four distinct canonical types of coordinates 
(not specific coordinates) that can be used to represent 
Si . In what follows we write out in each of these cases. 
We do not write out quantities that follow algorithmically 
from u°' . 



A. Kruskal-Szekeres {a = c = Q) 



The distinguishing characteristic here is 
A = 0. 

We can always choose 6 < so that 

1 



2hu 

and so it follows from © that 



T >0, 



G 



ABb^' 



(8) 
(9) 

(10) 



B. Diagonal (6 = 0) 

Here a < 0, c > and 

Be = Ga. 

With 

V = A^a + 4:B^c 
it follows that Da > 0, 



2 a ^ V Va 



and 



Va 



(11) 
(12) 

(13) 
(14) 



The requirement of "comoving" coordinates is = <^ 
B = G = 0. 



C. Bondi (c = or a = 0) 

It is sufficient to consider the case c = 0. Now 6^0 
and for 6 > a;^ is an "advanced" time and for 6 < a;^ 
is a "retarded" time. Now 



Ab = Ga. 

With 

£ = Aa'^+ iBab 
it follows that sign{£) — sign{A), 

1 4, A 



and 



(15) 



(16) 



(17) 



(18) 



The condition of "comoving" coordinates v? = requires 
a < and B = 0. Conversely, if i? = then v? = for 
a < and = —^fajb for a > 0. 



With 



D. a / 0, 6 / and c / 



7^ = A^c - 2A6C + 2CSc + 2(7^0 (19) 



and 



S = 2AbcB-2AbGa+B'^c^-ABGb^+2BcGa+G^a^+caA^, 

(20) 

it follows that 



12 _ -7^ + sign(b)y/n'^ - ASG"^ 
' " 25 



(21) 



and 



^2 ^ -bu^ + stgnib)^ibu^r^cia{u^r + l) ^^2) 



The condition of "comoving" coordinates = requires 
the same conditions as in the previous case. 



IV. PHENOMENOLOGY 

The discussion in this section requires no specification 
of coordinates on Si. 
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A. Decomposition 

Whereas condition ||SJ) determines the flow on Si, it 
is of interest to Icnow when the flow reduces to that of 
a perfect fluid (isotropic pressure (including bulk stress) 
and zero shear stress), or conversely, to be able to state 
that a given spacetime can not be compatible with a per- 
fect fluid source. The familiar case is entirely obvious: If 
the Einstein tensor is diagonal with three equal compo- 
nents, then the metric is consistent with a perfect fluid 
in comoving coordinates. Such a circumstance is, how- 
ever, a property of the spacetime and the coordinates in 
which it is exhibited, more general cases are explored in 
this section. With u" known, it is possible to reconstruct 
phenomenological parameters associated with a decom- 
position of the energy-momentum tensor directly from 
the metric. This decomposition is not always unique. 

A phenomenological fluid interpretation of Q , by way 
of Einstein's equations, follows from the well known 
Eckart relation |lC| which gives 

{VaT+ Tu^VpUc) ^0 (23) 

where k is the thermal conductivity and T the temper- 
ature proflle. The flows considered here can therefore 
be considered non-conducting {k = 0) setting aside con- 
trived functions T. The energy-momentum tensor is de- 
composed here in the form 11] 

= pu"up+pin°'np+p2Sp+p2{u''up-n°'n/3)-2riap, 

(24) 

where cr^ is the shear associated with u" and rj the phe- 
nomenological shear viscosity. It follows from ^ that 
G3 — G\ for the spacetimes considered in this paper. 
With Gl = G| it follows from ^ that either erf = cr| 
or we must set 77 = 0. Note that H24|l distinguishes 
the shear stress from an anisotropic pressure. These 
are sometimes combined. For example, in the comov- 
ing frame in a spherically symmetric spacetime, it follows 
that r| = diag{p, PI, P2, P2) where PI = and 
P2 = P2 — 2?7(T3. Such a combination is not, in general, 
possible outside the comoving frame. Since anisotropic 
pressures do not arise solely due to shear stresses {e.g. 
in a static spherically symmetric spacetime, ct^ = but 
Pi 7^ P2 in general), the full decomposition is used 
here. 



B. Systems of Equations 

We set up systems of equations to be solved simulta- 
neously for the functions (p,pi,p2,?7) in terms of scalars 
that follow algorithmically from the metric. It is not 
possible to build more than three independent equations 
in an attempt to solve for (p, pi,p2, r]) since, with ||SJ , 
there are only three independent scalars that can be con- 
structed from the set {G^,u" ,na) |3l- 



1. General Spacetimes 

To proceed in a manifestly invariant way we construct 
scalars from the set {G^,u" ,na) that are linear in G^. 
These are: G^u'^np (used in condition Q), G'^u°'u/3, 
G^n'^np and G^ = G. With ^ it follows that 

G^8n{-p + pi + 2p2), (25) 

and 

G^u^up EE Gl = Snp. (26) 
In all cases we take (|26|l as the definition of p. Further, 

Gin^Ufj = G2 = 87r(pi - 277A) (27) 

where 

A = o-f n"n/3. (28) 

Rearrangement of H25|) . H2t)|) and (|27|l gives 

487r?7A = G + Gl - 3G2 + 167rP (29) 

where P = pi —p2. If = or 77 = or P = then H25|l . 
(|26|l and 127|) form a complete set of equations. More gen- 
erally, however, if we attempt to solve for the complete 
set of parameters (p, pi,P2, v) four equations are needed 
[1^ . For the class of spacetimes considered here, it was 
shown in that higher-order invariants are not indepen- 
dent from the linear ones. Therefore, use of higher-order 
invariants will not break the degeneracy but merely gen- 
erate new syzygies (algebraic identities amongst invari- 
antly defined quantities). 



2. Restricted Spacetimes 

From 124|) and Einstein's equations it follows that for 
the spacetimes {Ai , g) considered here 

GlGl - {Gl - Gl){Gl - GX) ^G5 = {&7:)'P{p + P2), 

(30) 

where G'^ = G^ + IGTrrja'^. Although ^ holds in every 
{A4, g) without specific coordinates specified on Si, it is 
not manifestly invariant. Use of H3Ufl merely generates 
further (restricted) syzygies. 



C. Linear cases 

1. erf = 

If tj^ = and P+P2 7^ then equation with P ^ 
gives the Walker's pressure isotropy condition 13] 

GlGl^{Gl-Gl){Gl-Gl) (31) 
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which is here a necessary and sufficient condition for a 
perfect fluid. A manifestly invariant condition follows 
from (j23), lEH) and |(2ZJ which give 



where rj is arbitrary. If we set 77 = then the condition 
(|34|l is a necessary and sufhcient condition for a perfect 
fluid. For other choices of r/ the fluid is imperfect. 



Pi 



91 

8^' 



and 



Clearly 



P2 



G + GI-G2 
16^;^ ■ 



G + Gl = 3G2 



(32) 



(33) 



(34) 



is a necessary and sufficient condition for a perfect fluid 
including the exceptional case G + 3G1 = G2 not covered 
by . Equation is not sensitive to the presence of a 
cosmological constant term in the Einstein field equations 
and plays a central role in what follows 0| 

In all of what now follows up to section V we assume 



2. ri = 

The decomposition (|24|) is consistent with some space- 
times if and only if r/ = 0. Some examples of this are 
shown in section Ivl If 77 = then equations ()32|l . H33f) 
and H34|l hold as in the previous case. 



3. pi = p2 = p, A ^ 
Equations (j^SJ , if^ and lf77|l now give 



P 



G + Gl 
247r 



and 



rj = 



G + Gl - 3G2 
487rA 



(35) 



(36) 



so that is once again a necessary and sufficient con- 
dition for a perfect fluid. The case A = is equivalent 
to the case — 0. 



4. pi/p2, A/0 



Equations (^3)) , if^ and lf77|l now give 
G2 

pl = — + 277A, 



and 



„ G + G1-G2 



(37) 



(38) 



D. Non- uniqueness of the source 

In the last case above, it is not possible to solve for 
a unique set {p,pi,p2,ii) as only three invariants arc in- 
dependent for the type of spacetimes considered in this 
paper. Another way to see this directly is to observe that 
substitution of the expressions for p, pi and p2 as given 
by H2(j|) , H37|l and H38|l into the energy- momentum tensor 
(|^ and multiplication by Stt reproduces the Einstein 
tensor [l^. This of course all derives from the fact that 
in a canonical frame there are at most three independent 
components of the Einstein tensor for the spacetimes con- 
sidered. As we show in section Ivl the application of the 
framework to some given spacetimes known to represent 
perfect fluid solutions but where A 7^ shows that there 
are other imperfect fluid sources possible. This is always 
the case when the perfect fluid condition (|34|) holds and 
A ^ 0. For 77 = the fluid is perfect, for other choices 
of r] the fluid is imperfect. For example, the Lemaitre- 
Tolman-Bondi metric ("dust") is given by [l6|[^ 



-{dtf 



{R'{t,r))^{dr) 



1 + fir) 
along with the constraints 



R{t,rfdn^, (39) 



R{t,r) = J2 



m{r) 



fir), 



Rit,r) = - 



m(r) 
R^it,r)' 



(40) 



(41) 



and 



R {t,r) 



-1 , , m (r) m(r)R (t,r) 



2 TO {r)R{r,t) ~2m{r)R (r, t) + / (r)i?2(r,i) 



2i?(r, t) ^i2mir) + fir)R{r,t))R{r,t) 



where = -S- and 

or 



(43) 



Condition H34|) holds so the 



source is consistent with (but not necessarily) a perfect 
fluid, in fact simply dust (since G2 = = G + Gl) with 



TO (r) 



47ri?2(r, t) R'ir,t)' 



(44) 



The metric (|39|) (with the given constraints) is also con- 
sistent with an imperfect fluid withpi ~ 2riA, p2 = — ??A, 
77 arbitrary and 
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-2 m (r, t)R{r, t) + 6 m{r)R (r, i) - / (r)i?2 (r, t) + 2 R (r, t)f{r)R{r, t) 



3R' (r, t) )i?3/2 ^2m(r) 
I 



(45) 



as is easily verified. The fact that the Lemaitre metric 
need not he considered as dust is not a new result [l7l| . 
Here, the degenerate case (A 0) is shown to be generic. 



V. EXAMPLES 

We provide here examples in various coordinate types 
in order to illustrate the results obtained. 



A. Kruskal-Szekeres coordinates 

Aside from the Kruskal-Szekeres metric, little use has 
been made of double null coordinates. The example used 
here is the Einstein-de Sitter universe 113 



dsl^ = C^{u + v)\-dudv + ^" 



(46) 



where C is a constant and dVl^ is the metric of a unit 

:(1, 1,0,0), ai = 



sphere. It follows that — -pr, — i— 
and H34|l holds so that the fluid is necessarily perfect (in 
fact simply dust since G2 = = G + Gl). 



B. Bondi coordinates 



The Bondi metric [IE 



ds 



M 



c^iw, r)f{w, r){dwY ± 2c(w, r)dwdr 



r^dfl^ 
(47) 

in advanced (+) or retarded (— ) w has A — C~Q,B^O 
for ^ = and ^ ^ and so there is no non-conducting 
fluid source of ()47|l under these conditions. The Vaidya 
metric poj (corresponding to a null flux) provides a fa- 
miliar example. The metric j21| 

dsl^ = -2H{u, r){duf - 2dudr + ur^'^{{dxf -f {dyf), 

(48) 

where H{u, r) = (r/w+fcr"u(2-m)/(m-i))/2 > and n = 
m(m — l)/2 is necessarily comoving. Moreover, erf 7^ 
and (|34|l holds so the source is consistent with (but not 
necessarily) a perfect fluid. If we set ry = then the fluid 
is perfect but for other choices of 77 the fluid is imperfect. 



1. K=Q 

The Robertson- Walker metric 

a{tf{drf 



1 — fcr^ 



(49) 



gives erf = and l|34ll holds. It follows that the fluid is 
necessarily perfect 22]. Equations H26() and H27() repro- 
duce Friedmann's equations. In contrast, the Kantowski- 
Sachs metric 1231 



-{dtf 



a(tf{dr) 
1 — kr'^ 



dVL^ 



(50) 



gives erf = but (|34|l never holds. It follows that the 
fluid can never be perfect. For the general spherical static 
metric 

dsl, = -e'^(^){dtf + - 



(51) 



2m,{r)/r 

whereas erf — 0, (|34ll does not in general hold. The fluid 
has anisotropic pressure, the perfect fluid being a special 
case 113. The metric 123] 

ds^M = -{dtf + R{tf{{drf + sin{rf{dzf + f{rf{dcj,f) 

(52) 

with the constraint 2RR + (7?)^ + 1 = Q where ■ = ^ 
has erf = = 0. The condition G^ = gives /(r) = 
cos{r + A) where A is constant. Condition 134|) does not 
in general hold. Rather, P2 — ^ but pi = only for 
yl = 0. 



2. rj = 



The metric m 



ds]^ = S{t)-^'^C{x)-^-^-^{-{dtf + {dxf) 

+S{t)C{xr{T{tr{dyf + T{t)--{dzf){hi) 

has G% — — Gl and so is obviously a perfect fluid 
in comoving coordinates. However, er| ^ erf and so the 
decomposition H24() holds only for 77 = 0. Condition H34|l 
holds in agreement with the obvious. 



D. Non-comoving diagonal coordinates 



C. Comoving diagonal coordinates 

The next examples are in diagonal coordinates and 
have i? = G = so that is necessarily comoving. 



Few examples are available in non-comoving coordi- 
nates. From the pioneering work of McVittie and Wilt- 
shire we note for example, that their solution (6.12) 



ds 



M 
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where z = e(^^ — i^)/Coj ^-^id /^(-s^) is an undetermined 
function of z with /J^^ — /J^ ^ g, has = and the 
condition holds so the source is necessarily a perfect 
fluid. Their solution (7.20) 

dsXi = expiAe//^ + Bz/l - 2eLt){~{dtf + (duj)'^ + dfl^) 

(55) 

where z = lo + et, e = ±f and A, _B and L are constants, 
has (T^ 7^ 0. Condition 1)34(1 holds so the source is con- 
sistent with a perfect fluid (if 77 = 0). Their solutions 
(6.21) and (8.11) also have shear and (|34l) holds so the 
solutions are compatible with a perfect fluid source. As 
with the Davidson metric (|48|l . McVittie and Wiltshire 
solutions (6.21), (7.20) and (8.11) are also consistent with 
an imperfect fluid with rj ^ 0. 

From the more recent work of Senovilla and Vera 0, 
for example, their solution (40) 

has = Gl but 7^ erf and so the metric (|56|) is con- 
sistent with the decomposition H24() only for r/ = (the 
same holds for their solutions (38) and (41)). Equation 
(|34|l holds and so the metric (|56|l necessarily represents 
a perfect fluid. 

VI. SUMMARY 



is a necessary and sufficient condition for a perfect 
fluid. 

• iii. If A 7^ then pi is given by 1(37(1 and p2 is 
given by ijSHl) and ?7 is a freely specified function. 
This is a generic degenerate case. If condition ((34(1 
holds then the fluid is compatible with a perfect 
fluid source. If 77 = the fluid is perfect. For other 
choices of rj the fluid is imperfect. 

Furthermore, the derived covariant perfect fluid condi- 
tion can be used to study exact solutions of Einstein's 
equations as well as deriving new families of solutions 
|27j . Examples, in various coordinates, including space- 
times with and without unique interpretations have been 
provided. The procedure developed has been imple- 
mented in a computer algebra program described in the 
appendix. 
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It has been shown that a flow (timelike congruence u") 
in any type Bi warped product spacetime is uniquely and 
algorithmically determined by the condition of zero flux 
{G^u'^np = 0). Explicit forms of have been written 
out for canonical representations. With u" known, a phe- 
nomenological interpretation of the spacetime {Ai , g) in 
terms of a non-conducting fluid follows. The following 
cases are delineated 

• i. If 77A = then the condition ((34(1 is a necessary 
and sufficient condition for a perfect fluid. 

• ii. If pi = P2 = P and A 7^ then p is given by 

and 77 is given by The condition 



APPENDIX A: GRSOURCE, A COMPUTER 
ALGEBRA PROGRAM 

GRSource runs under the system GRTensorll 
which in turn runs under the system Maple. The program 
is called with a spacetime metric name as an input. GR- 
Source calculates the velocity field and related physical 
quantities. Moreover, it performs an algorithmic analysis 
on the nature of the fluid. The results of the evaluated 
quantities are displayed, followed by a summary report 
on the analysis of the possible fluid sources for the input 
spacetime. An example follows: 



> restart; We start a GRTensorll and GRSource session 



> grtw(); 



GRTensorll Version 1.80-pre2 (R6) 
Developed by Peter Musgrave, Denis Pollney and Kayll Lake 
Gopyright 1994-2003 by the authors. 
Latest version available from: http:/ /grtensor. or^ 



> qload(kantosachs); 



GRSource Package Version 1.00 
Developed by Mustapha Ishak and Kayll Lake, (c) 2002-2003 
Help available via ?grsource 
Usage: Load a metric and enter the command source(metricname) ; 



Default spacetime = kantosachs 



For the kantosachs spacetime : 
Coordinates : x(tip) 
x'' = [r, 0, 0, t] 
Signature = 2 
Line element : ds^ 



^^2^ &{t)^ dr^ ^ , +a(i)2sin(^)2 d<t>'' - dt' 

1 — A H 



Kantowski — Sachs metric, J. Math. Phys. 7, 443 {1966) 



> source(kantosachs); 



A general velocity field will be generated automatically from the metric. 
Please answer the following or enter exit anytime to stop the session. 

Enter the timelike or null coordinate, for example: t; > t; 
Enter the spacelike or null coordinate, for example: r; > r; 
Enter the spacelike coordinate, for example: theta; > theta; 
Enter the spacelike coordinate, for example: phi; > phi; 

For the kantosachs spacetime : 



u" = [0, 0, 0, 1] 
uSq = -1 
flux = All components are zero 
ShearTensor"' i, = All components are zero 

Expansion scalar : expsc 



Thetalu] = 3 ^^^V 



Delta = All components are zero 
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8 a(t)2 7r 



a(t)2 TT 



1 (Aa(t))2 + 2aft)(^a(t)) 
8 a(t)2 7r 

1 



PFCondition = 2 



a(t)2 
Report : 

eta * Delta = but PFCondition was not simplified to zero, 

the fluid has anisotropic pressure. 

Further simplifications can he applied to the objects calculated 

using the commands gralterQ and grdisplay{). 

I 
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